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Abstract 



The aim of this paper is to study a generalization of the two colors multitype contact process 
intended to mimic an example of interspecific competition called allelopathy. Roughly speak- 
ing, the strategy of the first species is to inhibit the growth of the second one by freezing 
the sites it has colonized. Our main objective is to prove that this spatial model exhibits 
some phase transitions comparable with those of the classical multitype contact process. In 
particular, we show that depending on its birth rate each species may conquer the other one. 

1 Introduction 

Allelopathy denotes a process involving secondary metabolites produced by plants, micro-orga- 
nisms, viruses and fungi that influence development of biological systems. Typically, the strategy 
of the so-called inhibitory species involved in such a process consists in increasing its compctitivity 
by inhibiting the growth of rival species, called susceptible species. The phenomenon is now well 
known and was observed in many ecosystems. Oueslati (2003) showed for instance that leaf 
extracts of Triticum durum L. (durum wheat) depresses the germination rate and radicel length 
of Hordeum vulgare L. (a barley variety) and Triticum aestivum L. (a bread wheat variety). 
Leaves of Pueraria thunbergiana also possess allelopathic activity by secreting a substance called 
xanthoxin that inhibits the root growth of Lepidium sativum L. (cress) seedlings. See Kato- 
Noguchi (2003). Allelopathy cases were also frequently observed in aquatic environment between 
marine algal species such as Alexandrium that has inhibitory effects on growth of Gymnodinium 
mikimotoi, Scrippsiella trochoidea and Chaetoceros gracile. See Arzul et al. (1998). Usually, these 
phenomena are characterized by the formation of empty spaces due to the production of toxins 
that prevent the susceptible species from setting up. 

The model we introduce to mimic this phenomenon is a continuous time Markov process whose 
state at time t is a function £ t : Z d — > {0, 1, 2, 3}. A site x in Z d is said to be occupied by the 
inhibitory species or blue particle if ^(x) = 1, occupied by the susceptible species or red particle if 
£(x) = 2, and empty otherwise. In the third case, x will be called a free site if t;(x) = and a 
frozen site if ^(x) = 3. To describe the evolution rules, we now let Af C Z d be the set of y € Z d 
such that || 2/ II < f where r is a positive constant and || ■ || some norm on R d . In other respects, 
we denote by f. L (x, £) the fraction of neighbors of x occupied in the configuration £ by a particle 
of type i, where the neighbors of x refer to the translated set x + Af. With the notations we have 
just introduced, we can formulate the transition rates as follows. 

1 Ai/ifoO 1 > 3 1 

3 — > 1 Xifi(x,£) 3^0 7 

.— 2 \2h(x,£) 2^0 1 



In particular, the model can be seen as a generalization of the two colors multitype contact process 
(see Neuhauser (1992)) in which the strategy of the blue particles is to inhibit the spread of the 
red ones by freezing the sites they have colonized for an exponentially distributed amount of time. 
Reciprocally, the multitype contact process can be considered as the extreme case 7 = 00, i.e. the 
transition from 3 to is instantaneous. 

We now formulate our results and construct step by step the phase diagram of the process. First 
of all, if we start the evolution from the measure 82 that concentrates on the all 2's configuration, 
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the process is called the basic contact process with parameter A 2 . In such a case, we know that 
there exists a critical value A c G (0, oo) such that the process converges in distribution to the all 
empty state if A 2 < A c , and to a measure /i2 that concentrates on configurations with infinitely 
many 2's otherwise. See, e.g., Liggett (1999). If we start from the all l's configuration, we have the 
same result. That is, if Ai < A c then £ t So while if Ai > A c then £ t converges to the stationary 
measure v\ = lim t ->oo#i S t , where St denotes the semigroup of the process. Moreover, since 3's 
do not disturb l's, it is not difficult to prove that u\{^{x) = 1) = /xi(£(x) = 1), where fix is the 
upper invariant measure of the basic contact process. 




Fig. 1 Configurations of the nearest neighbor multitype contact process with frozen state at time 
t — 50 on the square {0, 1, • • ■ , 199} 2 . Left picture: \i — M — 1,96 and 7 = 0,05. Right picture: 
Ai = 1,96, A2 = 2,88 and 7 = 0,05. The black points refer to the inhibitory species, the pale gray 
ones to the susceptible species, and the dark gray ones to frozen sites. 

To study the competition between both species, we now start the evolution with infinitely many 
l's and 2's and suppose that Ai > A c and A2 > A c . We first set 71 < 72 and denote by £ t ' 
the allelopathy model with parameters Ai, A2 and 7^. Then, if we think of the processes as being 
generated by Harris' graphical representation, we may run £ t x and £ t 2 simultaneously starting from 
the same initial configuration in such a way that £ t has more l's and fewer 2's than £ t , i.e. for 
any x G Z d , if (a;) = 2 then £ t 2 (a;) = 2, and if £ t 2 (a;) = 1 then (x) = 1. The same argument 
implies that the process is also monotonous with respect to each of the parameters Ai and A 2 . 

Theorem 1 We set t * = { x € Z d ; £t{x) = i } the sites occupied at time t by a type i particle. 
Then, for i = 1,2, the survival probabilities P (9 t * 7^ for all t > 0) are monotonous with respect 
to each of the parameters Ai, A2 and 7. 

In particular, by comparison with the multitype contact process (7 = 00), it looks clear that the 
l's win as soon as Ai > A 2 , i.e. £ t => v\. See Theorem 1 in Neuhauser (1992). More precisely, we 
have the following result. 

Theorem 2 We assume that £0 is translation invariant. If \\ > A 2 or Xi = X 2 and d > 3, then 
£i => v\ the limit starting from the all 1 's configuration. 

The proof of Theorem 2 for Ai > A2 simply relies on a coupling argument. On the other hand, 
if we now assume that Ai = A2 and d > 3, comparison with the multitype contact process just 
implies survival of l's for any 7 > 0. See Neuhauser (1992), Theorem 3. In this case, the key of 
the proof is duality. The ancestry of a given point (x, t) in Z d x R+ exhibits a tree structure that 
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can be divided, as for the ecological succession model (see Lanchier (2003)), into several layers, 
due to the transition 1 — > 3. The idea is that the second layer provides enough l's to freeze the 
path of the first ancestor at infinitely many points, which blocks 2's from determining the color of 
(x,t). Our approach is roughly the same as in Lanchier (2003), Sect. 4 and Sect. 5, and we think 
that it may be extended to a larger class of spatial models. 

If we now focus on the case Ai < A 2 , it is not clear that the 2's win since the dynamics is 
favorable to the l's. Nevertheless, by rcfcring again to the extreme case 7 = 00, we can prove that 
for fixed Ai < A 2 the 2's win in d = 2 and for 7 sufficiently large. 

Theorem 3 We suppose that d = 2 and \\ < A 2 . Then, there exists a critical value j c e (0, 00) 
such that £ t \i 2 for any 7 > -f c . 

The proof of Theorem 3 is quite simple. It essentially relies on the rescaling argument introduced 
in Durrett and Neuhauser (1997), Sect. 3. The idea is to prove that the argument resists to 
small perturbations of the process. That is, if 7 > -f c then, with probability close to 1, the al- 
lelopathy model behaves like the multitype contact process inside a given bounded space-time box. 

To draw the phase diagram of the process, we now define the critical value A 2 (7, Ai) as the 
infimum of A 2 > such that the l's die out, with the convention inf = 00. A straightforward 
application of our theorems then implies that A 2 (7, Ai) J, Ai as 7 | 00 and A 2 (7, • ) is increasing 
on R+. In conclusion, the picture we finally obtain is given by Fig. 2. Unfortunately, we don't 
know if the 2's may coexist with the l's for Ai < A 2 < A 2 (7, Ai) (gray part of the picture), and if 
they can win for any Ai > A c , that is if A 2 (7, Ai) < 00 for any Ai > A c . To have an idea of the 
answers for these questions, we can resort to the mean field theory. More precisely, if we pretend 
that the functions Ui(x,t) — P (£t(x) = i) do not depend on the site x G Z d , the evolution of the 
process can be formulated in terms of the following ordinary differential equations. 

iio = u 2 + 7113 - XiUqUi - X 2 U U 2 
u[ = Ai uq ui + Ai uq U3 — ui 
u' 2 = A 2 uo u 2 - u 2 

u' 3 = Ml — Al Ul U3 — 7U3 

Let n — { u ; Ui > 0, uq + u\ + u 2 + 113 = 1} be the collection of values we are interested in. 
For fixed 7 > 0, we also denote by lo\ (resp. lo 2 ) the set of parameters Ai and A 2 such that Ai > 1 
and 7 A 2 < Ai (Ai + 7 — 1) (resp. A 2 > 1 and A 2 > Ai). First of all, a straightforward calculation 
shows that the ODE has a nontrivial fixed point u on the boundary u 2 — if and only if Ai > 1, 
where nontrivial means ^ (1,0,0,0). Moreover, by studying the eigenvalues of the linearization 
at point u of the ODE, we can prove that the equilibrium u is attracting if (Ai,A 2 ) G u>i, and 
unstable otherwise, that is the linearization has an unstable direction that points into int(£l). In 
the same way, if A 2 > 1 there exists a nontrivial equilibrium v on the boundary v\ = V3 — that 
is attracting if (Ai,A 2 ) £ u> 2 , and unstable otherwise. Finally, one can prove that the ODE has 
a fixed point belonging to int(fi) if and only if (Ai, A 2 ) 6 uji n uj 2 . In other words, coexistence is 
possible for the mean field model as soon as Ai < A 2 < 7 _1 Ai (Ai + 7 — 1). 

Typically, if the range of the interaction r or the dimension d are sufficiently large, one expects 
that the process has the same properties as the mean field model. In particular, relying on the 
instability of u for (Ai, A 2 ) 6ui, we believe that for any Ai > A c the 2's win for A 2 large. However, 
in view of the properties of the mutitype contact process (See Neuhauser (1992)), we don't think 
that coexistence may occur in an open set of the parameters Ai and A 2 . In conclusion, we formulate 
the following conjecture. 

Conjecture 4 // the range of the interaction r > ro or the dimension d > do, then for any Ai 
there exists A 2 (7, Ai) < 00 such that the 1 's win if A 2 < A 2 (7, Ai) and the 2's win if A 2 > A 2 (7, Ai). 
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Fig. 2 Phase diagram. 



2 Graphical representation and duality 

We begin by constructing the process from a collection of Poisson processes in the case Ai > A2 . 
For x, y £ Z d , x - y e N, let {T£'» ; n > 1}, ; n > 1} and {V* ; n > 1} be the arrival times 
of Poisson processes with rates Ai/ card A/", 1 and 7 respectively. At times T£ ,y , we draw an arrow 
from x to y, toss a coin with success probability (Ai — A 2 )/Ai, and, if there is a success, label the 
arrow with a 1. If x is occupied by a red particle, that the arrow is unlabelled and that y is free, 
the particle will give birth through this arrow. If x is occupied by a blue particle and that y is 
free or frozen, the site y will be painted in blue. Now, at times U%, we put a cross x at x. The 
effect of a x is to kill both species, i.e. a red particle gives way to a free site and a blue one to a 
frozen site. Finally, at times V*, we put a dot • at x to indicate that a frozen site becomes free. 
A result of Harris (1972) implies that such a graphical representation can be used to construct the 
process starting from any initial configuration £0 S {0, 1, 2, 3} z . For an example of realization of 
the process, see Fig. 3. 

After constructing the graphical representation, we can now define the dual process. We say 
that two points (x,s) and (y,t) in Z d x R + are connected or that there is a path from (x,s) to 
(y, t) if there exists a sequence of times so = s < si < s 2 < • • • < s n < s„+i = t and spatial 
locations xq = x, x±, X2, ■ ■ ■ x n = y so that 

(1) for 1 < i < n, there is an arrow from to Xi at time Sj and 

(2) for < i < n, the vertical segments {x^ x (sj, s i+1 ) do not contain any x's. 

If (1) only is satisfied, we say that there is a weak path from (a;, s) to (y, t). Finally, we say that 
there exists a dual path from (x,t) to (y, t — s), < s < t, if there is a path from (y, t — s) to 
(x,t). In other words, dual paths move against the direction of time and arrows. We then define 
the dual process by setting 




{ y G Z d ; there is a dual path from (x, t) to (y, t — s)} 
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Fig. 3 Harris ' graphical representation. The black lines refer to blue particles, the pale gray 
ones to red particles, the dark gray ones to frozen sites, and the dotted ones to free sites. 

for any < s < t. First of all, wc can observe that {(£ s <!C,t) , s) ; < s < t }, which is the set of 
points in Z d x [0,t] that are connected with (x,t), exhibits a tree structure. As for the mutitype 
contact process, such a structure allows to equip the dual process £ s <x,t> with an ordered relation 
for which the members are arranged according to the order they determine the color of (x, t). See 
e.g., Neuhauser (1992), Sect. 2. From now on, the tree 

T = {(£"''\a);0<a<t} 

will be called the upper tree starting at (x,t) and the elements of £ s x,t the upper ancestors. We 
will denote by Ci^'' (n) the n-th member of the ordered ancestor set. For an example of ancestor 
hierarchy, see Fig. 4. If there is a weak path from (y, r) to (x, t), the tree starting at (y, r) will be 
called a lower tree and the elements of £ s <!/ ' r) the lower ancestors. Finally, the first upper ancestor, 
that is £ s (x,t> (1)) wm be called the distinguished particle. 

To conclude this section, we now describe in greater detail a general method to determine the 
color of (x,t) in the case Ai > A2. Contrary to the multitype contact process, the state of some 
sites (free or frozen) depends on the structure of the lower trees. So, for more convenience, we first 
assume that this information is known and say that an arrow is forbidden for the red if its target 
site is frozen or if it is a 1-arrow. We will see further that our results can be proved by focusing 
only on the first two layers of the tree structure. First of all, if the distinguished particle lands 
on a blue site, it will paint (x, t) in blue. If it lands on a red site and does not cross any arrow 
forbidden for the red, it will paint (x, t) in red. On the other hand, if it lands on a red site and 
crosses at least one arrow forbidden for the red, we follow the path it takes to reach (x, t) until we 
first meet a forbidden arrow. Then, we look backwards in time starting from the location where 
this arrow is attached and discard all the ancestors of this point, that is the next few members of 
the ancestor hierarchy. Finally, we start over again with the next remaining ancestor, and so on. 
In the same way, if the distinguished particle lands on a free or frozen site, we start afresh with 
the second ancestor. 

For instance, in Fig. 4 the distinguished particle lands on a red site and crosses an arrow that 
points at the frozen site x — 3 so fails in painting (x, t) any color. The same holds for the second 
ancestor that takes a 1-arrow to reach (x,t). In other respects, the third ancestor lives at time 
on the tree starting at the location where this 1-arrow is attached, so we look at the fourth 
ancestor. This last one finally succeeds in painting (x, t) in blue. 
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3 Proof of Theorem 2 

In this section, we will prove Theorem 2 beginning with the case Ai > A 2 . The strategy of 
the proof is quite simple. It relies on a basic coupling between the allelopathy model and the 
multitype contact process. First of all, recall that the multitype contact process can be considered 
as the extreme case 7 = 00, i.e. the transition from 3 to is instantaneous. In that case, it is 
known that if £0 S {0,1, 2} z is translation invariant and Ai > A2 then starting from infinitely 
many l's and 2's the first species win the competition, i.e. £ t => \i\ the upper invariant measure of 
the basic contact process. See Theorem 1 in Neuhauser (1992). Now the intuitive idea is that the 
smaller 7 is the more the blue particles are competitive. To specify this, we couple both processes 
together considering two identical copies G\ and G2 of the graphical representation. In the first 
copy, the particles evolve according to the rules given above. In the second one, the x 's permute 
both species into free site, whereas the dots have no more effect on the particles so G2 can be 
seen as the graphical representation of the multitype contact process. Then, we may run both 
processes and £ t 2 simultaneously on G\ and G 2 starting from the same configuration in such a 
way that has more l's and fewer 2's than £ t 2 , i.e. for any point (x, t) in Z d x R + , if £ t * (x) = 2 
then £ t 2 (x) — 2, and if £ t 2 (x) = 1 then (x) = 1. This together with Theorem 1 of Neuhauser 
(1992) implies the first part of Theorem 2. 
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The extension of this result for Ai = A2 and d > 3 is more difficult. In this case, the strategy 
consists in proving that if the upper tree T lives forever, the distinguished particle jumps infinitely 
often on a frozen site. To do this, we will focus on the second layer of the tree structure to show 
that, conditioned on survival of the upper tree, the number of blue particles that freeze the sites 
visited by the first ancestor can be made arbitrarily large. We will finally conclude by exhibiting 
an upper ancestor that will paint (x, t) in blue. 

From now on, we denote by A the common value of Ai and A2 and suppose that V lives 
forever. We can observe that such an event obviously occurs with positive probability since A is 
supercritical. For more convenience, we reverse the arrows and time by letting s = t — s. The 
main objective is to prove that the number of frozen sites visited by the first ancestor tends to 
infinity as t | 00. To specify this result, we start by defining more precisely our framework. First 
of all, we follow the path of the distinguished particle starting from (x, 0) and, for any n > 1, 
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denote by (3 n = (3 n (x,t) the n-th crossed arrow. See Fig. 5 for a picture. Now let z n = z n (x,t) 
and s n — s n (x,t) be respectively the starting site and the temporal location of the arrow (3 n (x,t) 
and denote by iV(a;, t) the number of arrows (3 n (x, t) that start at a frozen site, i.e. 

N{x,t) = card{n>l; = 3}. 

By construction of the sequence j3 n (x,t), it is clear that N(x,t) also denotes the number of frozen 
sites visited by the distinguished particle. In particular, the main result we have to prove can be 
formulated as follows. 

Proposition 3.1 If Xi = A 2 and d>3 then lim t -> oo N(x,t) = oo a.s. 

The intuitive idea of the proof is that the second layer of the tree structure provides enough 1 's 
to freeze the path of the distinguished particle at infinitely many points so that we can focus only 
on this one. To specify this point, we proceed in two steps. 

First of all, we denote by a n = a n (x, t) the arrival time of the first x located under the starting 
point of (3 n (x, t), i.e. 

& n (x,t) = min{[/^" such that U£ n > s n (x,t)} 

k>l 

and let T n — T n (x,t) be the lower tree starting at (z n ,a n ), that is 

r„ = {(y, s) e Z d x [<r„,oo) ; there is a path from (z n ,a n ) to (y,s)}. 

See Fig. 5 for a picture. We say that the n-th lower tree is favorable to the 1 's if the following 
two conditions are satisfied. 

(1) T n (x,t) lives forever and 

(2) the vertical segment {z n } x (s n ,a n ) does not contain any dots. 

As we will see further, the properties (1) and (2) will give us a good opportunity to freeze the site 
z n at time s n so the first step of the proof is to show that there exist infinitely many lower trees 
r„ that are favorable to the l's. 

Lemma 3.2 Let G n be the event that the lower tree T n is favorable to the l's. If the upper tree 
lives forever then P (limsup,^^ G n ) = 1. 

Proof. To begin with, denote by A n the event that T n lives forever and by B n the event that 
{z n } x (s„,<7„) does not contain any dots. The first step is to prove that for any n > 1, there 
exists a.s. an integer m > n such that A m occurs. To do this, we set T ni = T n and while T Uk 
is bounded we denote by r„ fe+1 the first lower tree that borns after T nk dies. Note that if A nk 
does not occur then r„ fe+1 is well defined and the event A nk+1 is determined by parts of the graph 
that are after T nk dies so that A nk and A nk+1 are independent. More generally, since the trees 
Tni, r„ 2 , • • • , r„ fc+1 are disjoint, the events A ni , A„ 2 , • • • , A Uk+1 are independent. Moreover, the 
probability that A nk occurs is given by the survival probability p\ of the basic contact process 
with parameter A starting from one infected site so that 

p(A c n n A r c l+1 n •••) < lim p ( A c ni n A c n2 n • • • n A c n ) 

K—>OC 
OO 

< U P ( A rJ = l™ (!-PA) fc = 
k=l 

as soon as A > A c . In particular, 

P(limsup,4„) = lim P(A n U A n+1 U •••) = 1. 
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This proves that, with probability 1, there exist infinitely many lower trees that live forever. In 
other respects, since a n — s n is exponentially distributed with parameter A, we can state that for 
any n > 1 

P(B n ) = P(a n -S n <V 1 z ") = A 7 - 1 (A + 7 )- 1 > 0. 
By independence, we can finally conclude that P ( lim sup^^ A n n B n ) = 1. □ 

To complete the proof of Proposition 3.1, we now consider for any n > 1 and s > a n the time 
translation dual process 

i^"' 5 "' = {yeZ d ; there is a path from (z n ,a n ) to (y,s)} 

and denote by j s (n) the associated distinguished particle, that is the first ancestor of (z n ,a n ). 
First of all, we can observe that if the lower tree T n lives forever then 7s (n) is well defined for 
any s > a n . Moreover, if we suppose that T n is favorable to the l's and that 7s (h) lands at time 
t on a blue site it is clear, in view of the condition (2) above, that z n will be frozen at time s n . 
In particular, letting T nk be a subsequence of favorable trees given by Lemma 3.2, the proof of 
Proposition 3.1 can be completed with the following lemma. 

Lemma 3.3 Let Cl s = { 7s (nfc) ; k > 1} and be the set of sites occupied at time s by a 1. If £ 
is translation invariant and d > 3 then, starting from infinitely many 1 's, 

lim card(f2 t fl a) = oo a.s. 

Proof. The key of the proof is transience. More precisely, a straightforward application of Lemma 
5.4 of Lanchier (2003) implies that lim^oo card(f2 t ) = oo a.s. as soon as d > 3. We then con- 
clude with Lemma 9.14 of Harris (1976) and translation invariance of the initial configuration. □ 

To conclude the proof of Theorem 2, we now restore the direction of time and arrows and 
construct a sequence of upper ancestors £ t (x '*' (k), k > 0, that are candidates for painting (x, t) in 
blue. The first member, that is Ct <x t> (0), is the distinguished particle £ t (x ' 4) (1). Next, we rename 
the sequence of frozen points (zk,Sk), k > 1, visited by the distinguished particle by going forward 
in time. For fixed k > 1, we now look backwards in time starting from the location where the arrow 
Pk(x, t) is attached and discard all the offspring of this particular point. We then define Q x,t (k) 
as the first upper ancestor that is left after discarding. We now set ( t — (k) ; k > 0}. By 

Proposition 3.1 and the fact that the upper tree T is linearly growing in time, the cardinality of 
C t can be made arbitrarily large by choosing t large enough. In particular, a new application of 
Lemma 9.14 in Harris (1976) gives us that 

lim P(( t -i ne; = 0) = 0. 

t — >oc 

where 6^ denotes the set of sites occupied at time s by a 1. Hence, there exists at least one 
candidate that lands on a blue site. Denote by Ct <x (ko) the first one in the hierarchy. Since by 
construction the arrow f3k {x, t) is forbidden for the red particles (we recall that at time Sk the 
site Zk Q is frozen), the upper ancestor £ t ( x ' l) (ko) will finally paint (zk al Sk a ) and so (x,t) in blue. 
This completes the proof of Theorem 2. 

4 Proof of Theorem 3 

In this section, we assume that d = 2, set Ai < A 2 and prove that there exists -f c > such 
that for any 7 > -f c the 2's win. In view of the dynamics, the survival of 2's for Ai < A 2 is not 
clear and tools as basic coupling and duality techniques fail in proving Theorem 3. We will first 
rely on the rescaling argument in Durrett and Neuhauser (1997), Sect. 3, applied to the multitype 
contact process, and then prove that taking 7 large enough does not affect too much the process. 
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We start by introducing the suitable space and time scales. We let L define a positive integer and, 
for any z — {z\,z 2 ) in Z 2 , set 

$(z) - (LzuLsb), B = [-L, L] 2 , B(z) = <5>{z) + B. 

Moreover, we tile B(z) with L"' 1 x L"' 1 squares by setting 

ttH = (L 0A Wl ,L ' 1 w 2 ), D = (-L°' 1 /2,L°' 1 /2} 2 , 

D(w) = tt(w)+D, I z = { w e Z 2 ; D(w) C B(z)}. 

We say that B(z) is good if B(z) is void of l's and has at least one particle of type 2 in each of 
the squares D(w) for w £ I z . For z = (z\, z 2 ) € Z 2 with Z\ and z 2 both even for even k, and Z\ 
and z 2 both odd for odd k, we say that (z, fc) is occupied if -B(z) is good at time kT, where T is an 
integer to be picked later on. Moreover, we require this event to occur for the process restricted 
to the space box [—ML, ML] 2 + $(z). We start by assuming that £ t is the multitype contact 
process, that is 7 = 00. 

PROPOSITION 4.1 (Durrett and Neuhauser) Assume that X 2 > \i and T = L 2 . For any 
5 > 0, there exist large enough L and M such that the set of occupied sites dominates the set of 
open sites in an M-dependent oriented percolation process with parameter 1 — 5. 

See Durrett and Neuhauser (1997), Proposition 3.1 and Lemma 3.7. To extend the comparison 
result to 7 > sufficiently large, we just need to prove that, with probability close to 1, the process 
behaves like the multitype contact process (i.e. none of the 2's is blocked by a frozen site) inside 
the space-time box 

$(z) + [-ML/ 3, ML/3} 2 x [0,T] 

See Lemma 3.7 in Durrett and Neuhauser (1997). Clearly, this occurs if we free all the target sites 
by setting a dot under each tip of arrow. So, letting t(x, t) be the number of arrows that point 
at x by time t and n{z) = $(z) + [—ML / 3, ML / 3] 2 , and decomposing according to whether 
l{x, T)>2 X 2 T or i(x, T) < 2 A 2 T, we obtain 

P ( any of the 2's is blocked ) < ^ P (t(x, T) > 2 \ 2 T) + 2 \ 2 T ^ P ( V? < Uf) 

< (2/3)ML(Ce- /3T + 2A 2 T( 7 ( 7 +1))- 1 ) < S 

for T and 7 sufficiently large and appropriate C < 00 and (3 > 0. Comparison with M-dependcnt 
oriented site percolation and Proposition 4.1 then imply Theorem 3. 

Acknoledgement. I would like to thank Claudio Landim, Olivier Benois and Roberto Fernandez 
for their advice, and Pierre Margerie for the biological interpretation. 
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